home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / science / sm32a.zip / LIBRARY / PLOT.LI < prev    next >
Text File  |  1994-08-11  |  699b  |  20 lines
  1. #    plot(y,x)
  2. #    plots a function of y=y(x)
  3. #    on xy-plane by default xmin:=-5, xmax:=5, ymin:=-5, ymax:=5
  4. #    e.g. plot(sin(x),x)
  5.  
  6. plot(y_, x_, xmin_, xmax_, ymin_, ymax_, color_) := block(numeric:=on,
  7.     dx:=(xmax-xmin)/getmaxx,
  8.     dy:=(ymax-ymin)/getmaxy,
  9.     xstep:=dx*10,
  10.     axis(xmin,xmax,ymin,ymax),
  11.     setcolor(color),
  12.     moveto(0,getmaxy-(subs(y,x=xmin)-ymin)/dy),
  13.     do( lineto((x-xmin)/dx,getmaxy-(y-ymin)/dy),x,xmin,xmax,xstep),
  14.     numeric:=off,
  15.     readchar,
  16.     local(dx,dy,xstep))
  17. plot(y_, x_, xmin_, xmax_, ymin_, ymax_) := plot(y,x,xmin,xmax,ymin,ymax,colorno(yellow))
  18. plot(y_, x_, xmin_, xmax_) := plot(y,x,xmin,xmax,-5,5,colorno(yellow))
  19. plot(y_, x_) := plot(y,x,-5,5,-5,5,colorno(yellow))
  20.